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Abstract 

The stochastic scenario of relaxation in the complex systems is presented. It is based 
on a general probabilistic formalism of limit theorems. The nonexponential relax- 
ation is shown to result from the asymptotic self-similar properties in the temporal 
behavior of such systems. This model provides a rigorous justification of the energy 
criterion introduced by Jonscher. The meaning of the parameters into the empirical 
response functions is clarified. This treatment sheds a fresh light on the nature of 
not only the dielectric relaxation but also mechanical, luminescent and radiochem- 
ical ones. In the case of the Cole-Cole response there exists a direct link between 
the notation of the fractional derivative (appearing in the fractional macroscopic 
equation often proposed) and the model. But the macroscopic response equations, 
relating to the Cole-Davidson and Havriliak-Negami relaxations, have a more gen- 
eral integro-differential form in comparison with the ordinary fractional one. 
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1 Introduction 



Experimental investigations surely have established the relaxation response 
of various complex systems (amorphous semiconductors and insulators, poly- 
mers, molecular solid solutions, glasses, etc.) to be non-exponential in nature 
p!p] . In particular, all types of the empirical functions used to fit the dielec- 
tric data exhibit the fractional-power dependence of the dielectric responses 
on frequency and time. It is worth noticing that this unique property is in- 
dependent on any special details of examined systems. In the past decades, a 
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considerable attention has been paid to find a theoretical explanation of the 
experimental results [3Pf5|6] . The main feature of all the dynamical processes 
in the complex systems is their stochastic background. In this framework one 
can expect that the macroscopic behavior of the complex systems is governed 
by "averaging principles" like the law of large numbers to be in force. How- 
ever, to develop the assumption is enough difficult. The point is that their 
macroscopic evolution is not attributed to any particular object taken from 
those forming the system. The problem of constructing an "averaged" object 
representing the entire relaxing system is not trivial. The description of the 
relationship between the local random characteristics of complex systems and 
the universal deterministic empirical laws is of great importance. 

Many-body effects play a vital part in such systems. No wonder that there ex- 
ists a direct relationship, suggested in literature |l|7|Sf^ . between anomalous 
relaxation and anomalous diffusion. One of more convenient languages for the 
description of anomalous diffusion is the continuous random walk (CTRW) 
theory. It occupies an important place for studying many physical phenom- 
ena. The notation of CTRW was first proposed by Montroll and Weiss in 1965 
|1U] . With their happy touch the CTRW generalized a simple random walk. 
Although the term "random walk" was introduced by Pierson in 1905, the 
formalism of simple random walks was known else in the XVII-th century. 
The random walk approach is based on the assumption that step changes are 
made through equal time intervals. This was a first approximation in vari- 
ous physical, chemical and economical models. For its turn the CTRW theory 
allows a random waiting time among subsequent random jumps. To sum up 
the long-term studies of this problem, the recent, mathematically excellent 
works of Meerschaert, Scheffier and Becker-Kern [11] have made its details ul- 
timate clear. On the one hand, this has allowed one to recognize the stochastic 
processes responsible for the anomalous behavior. On the other hand, the ap- 
proach proposes a description of the anomalous properties. In fact, it has es- 
tablished a close connection between the stable distributions (from the theory 
of probability) and the fractional calculus. This means that the nondifferen- 
tiable nature of microscopic dynamics of components in the complex systems 
can be transmitted to the macroscopic description of such systems in the form 
of fractional operators [T^P3] . Consequently, the CTRW method is very pop- 
ular for physical applications connected with anomalous diffusion, transport 
in disorder media, superslow relaxation, etc. (see the perfect reviews p31[T5] 
and references therein). 

In this paper we suggest the probabilistic approach to the analysis of evolu- 
tion processes. Our approach is based on the probabilistic formalism of limit 
theorems which provides tools to relate the local random characteristics of 
complex systems to the deterministic and universal relaxation laws regardless 
of the specific nature of the systems considered. We attempt to answer the 
following key questions related to the temporal evolution of complex systems: 
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o what does mean a self-similarity in the evolution of the complex stochastic 

systems (Section 2); 
o how to interpret theoretically the empirical deterministic relaxation laws 

(Section 3); 

o what characteristics of the internal structure of the complex systems stand 

behind the empirical responses (Section 4); 
o what connection is between the micro/meso/macrocopic dynamics of the 

relaxing systems and the macroscopic energy criterion (Section 5); 
o what role of the macroscopic response equations is for the description of the 

relaxation phenomena (Section 6). 

Finally, we discuss some alternative models. 



2 Self-Similarity of Complex Systems 

The simplest traditional interpretation of relaxation phenomena is based on 
the concept of a system of independent exponentially relaxing objects (for 
example, dipoles) with different (independent) relaxation rates [TB]. Since any 
macroscopic system consists of a finite number of objects, the approach gives 
a discrete set of relaxation rates. Without a doubt, the assumption may be 
valid for some relaxing systems. However, nether finite nor infinite exponential 
series with different real rates and contributions can result in the Havriliak- 
Negami (HN) and Kohlrausch- Williams- Watts (KWW) response laws exactly. 
The only opportunity to overcome this problem is a mathematical extension of 
this model. Really, the transition from the discrete "topology" of the relaxation 
rate set to the continuous one changes the situation in quality. By the integral 
transform (resembling the Laplace transform and replacing the exponential 
series) one can obtain the well-known empirical response functions. In this 
case the temporary relaxation distribution density is a continuous function 
and has no any narrow peaks that could be interpreted as a manifestation of 
separate subsystems (objects). 

The relaxation, following the Cole-Cole (CC) law, may be developed in a two- 
state system. Let N be the common number of dipoles in a dielectric system. 
If iVi is the number of dipoles in the state is the number of dipoles in 

the state J, so that = A^^ -|- A''^^. Assume that for t = the system is stated 
in order so that the states t dominate, namely 



where is the part of dipoles in the state ^, the part in the state J,. 
Denote the transition rates by w defined from microscopic properties of the 
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system (for instance, according to the given Hamiltonian of interaction and the 
Fermi's golden rule). In the case the kinetic equation describing the ordinary 
relaxation (Debye law) takes the form 

h^{t) -w{n^{t) -n^{t)} = 0, 

h^{t)-w{n^{t)-n^{t)} = 0, ^ ^ 

where, as usual, the dotted symbol means the first-order derivative. The steady 
state of the system corresponds to equilibrium with n-|^(oo) = nj^(oo) = 1/2. 
Clearly its response has an exponential character. However, this happens to be 
the case for such dipoles that relax irrespective of each other and of their en- 
vironment. If the dipoles interact with their environment, and the interaction 
is complex (or random), their behavior already will not be exponential. 

Assume that the interaction of dipoles with environment is taken into ac- 
count with the aid of the temporal subordination. Recall that in the theory 
of anomalous diffusion the notation of subordination occupies one of the most 
important places. So, a subordinated process Y{U(t)) is obtained by random- 
izing the time clock of a random process Y{t) by means of a random process 
U{t) called the directing process. The latter process is also often referred to as 
the randomized time or operational time [T7j. Generally speaking, the process 
Y may be both random and deterministic in nature. The anomalous diffusion 
theory studies, as a rule, the subordination of random processes. We intend 
to extend this approach to relaxation processes. 

Let the time variable be a sum of random temporal intervals Tj being non- 
negative independent and identically distributed so that the waiting times Tj 
belong to an a-stable distribution (0 < a < 1). Then their sum n~^^°'(Ti + 
T2 + ■ ■ ■ + Tn), n G N converges in distribution to a stable law with the same 
index a ^llj. To determine a walker position at the true time t, one needs 
to find the number of jumps up to time t. This discrete counting process is 
{Nt}t>o = max{n G N | J2i=i ^« — 0- Denote the continuous limit of {Nt}t>o 
by S{t). For a fixed time it represents the first passage of the stochastic time 
evolution above that time level. The random process is nondecreasing, and it 
can be chosen as a new time clock (stochastic time arrow) ^18j . The probability 
density of the process S{t) has the following Laplace image 

^^•^ Br 

where Br denotes the Bromwich path. This probability density has a simple 
physical interpretation. It determines the probability to be at the internal 
time (or so-called operational time) r on the real (physical) time t [18]. The 
function Fa{z) can be expanded as a Taylor series. Besides, it has the Fox' 



4 



H-function representation 




£fc!r(l-a(l + A;)) 



where r(a;) is the ordinary gamma function. In the theory of anomalous dif- 
fusion the random process S{t) is apphed for the subordination of Levy (or 
Gaussian) random processes pT|IT9] . The inverse Levy process S{t) accounts 
for the amount of time that a walker does not participate in the motion pro- 
cess [20j . If the walker participates all time in the motion process, the internal 
time and the physical (external) time would coincide. 



As was shown in [TH], the stochastic time arrow can be applied to the general 
kinetic equation. Then the equation describing a two-state system takes the 
following form 



■ D^n^it) - w {n^it) - n^it)} = 0, q < a < 1 



(3) 



where i)" is the a-order fractional derivative with respect to time. Here we 
use the Caputo derivative [2T|22] . namely 

1 } x^'^Ut) 

D°'x(t) = — / ^ — - — dr, n — 1 < a < n, 

^ ' r{n-a)J{t- r)"+i-" ' 



where x(")(t) = D"a;(t) means the n-derivative of x{t). The relaxation function 
for the two-state system is written as 

(pcc{t) = 1 - 2n^{t) = 2n^{t) - 1 = E^{-2wr), 



where Ea{z) = Y^'^=q z"' /T(l+na) is the one-parameter Mittag-Leffler function 
[23] . Feller conjectured and Pollard proved in 1948 that the Mittag-Leffler 
function Ea{—t) is completely monotonic for t>0, ifO<a<l. Moreover, 
Ea{—t) is an entire function of order 1/a for a > [T7]. It should be pointed 
out that the relaxation function under interest corresponds to the CC law. 

The analysis of this problem will not be complete, unless one consider it on 
the other hand. It turns out that the same result mentioned above can be also 
obtained by another way Pl24f25j . If the relaxation rate of the z-th dipole is 
equal to the value 6, then the probability, that this dipole did not change its 
initial orientation prior to an instant t, is 

Pr {Oi>t\Pi = h) = exp{-bt) for t > 0, b> 0. (4) 
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The random variable /3i denotes the relaxation rate of the i-th dipole and 
the variable 6i , the time needed for changing its initial orientation. Let 
and {6i} form sequences of nonnegative independent identically distributed 
random variables. Following |2S] and the law of total probability, one define 
the relaxation function 0i(t) for z-th dipole as a probability: 

oo 

Mt) = Pr {e^ >t) = J exp{-bt) dHp^h) , (5) 





where if/j. is the distribution function of each relaxation rate The form of 
a suitable function Hj^. should be found. In the system consisting of a large 
number of relaxing dipoles, the relaxation function (f){t) has to express in 
terms of the probability that the entire system will be without changing its 
initial state until t: 

(f){t) = lim Pr (.4^ min(0i, . . . , ^^v) > t), (6) 

N—>-oo 



where An is a. normalizing constant. Let us observe that the expression ([5]) is 
the Laplace transform of the distribution function if/3.(6): 

Pr {e,>t)=C{Hp^;t). 



Because of 9i being independent, we get 

Pr (n,inW,....9„) > ±-) = (p. > J-)) " = J-)) 



When N tends to infinity, the A^-th power of the Laplace transform of the non- 
degenerate distribution function Hjj- converges to a non-degenerate limiting 
transform, if and only if if/j. belongs to the domain of attraction of the Levy- 
stable law IT71I26II27I . 



In fact, the above limiting form is only determined by the behavior of the tail 
of Hjs.{b) for large b, i. e. by asymptotic properties of Hp^{b). The detailed 
knowledge of its other properties is not necessary. It is enough that the dis- 
tribution function Hp^ belongs to a domain of attraction of the Levy-stable 
law with the index of stability a. On the other words [T7], the necessary and 
sufficient condition for any x > is 

limi^-M=x-". (7) 
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This condition can be interpreted as a type of self-similarity. Really, for any 
X > and for large b 

Ft {/3i > xh) ^ Pr {/3i > b). (8) 

It is that the self-similarity is suggested as a fundamental feature of relaxation 
phenomena [7|28|29] . It should be stressed here that in the approach this 
conclusion arises from the pure probabilistic analysis, independently of the 
physical details of dipolar systems. Thus, it can be carried over to other similar 
cases of complex systems. 

The randomness of the relaxation rates (3i (1< i < N) is motivated by the 
fact that in the complex systems an object has not the only equilibrium state, 
but their states form a whole set of metastable substates. Their configuration 
changes in a very complicated way during their evolution. Each of the objects 
is locked into a substate, and the distribution of relaxation rates rejects any 
deterministic behavior of an individual object in the complex system. The 
total survival probability of the whole system has formally the same form 

oo 

Pr {Oi >t) = (exp(-At)) = J exp{-bt) dH^{b) (9) 



like However, now the form of d. f. is strictly fixed (i. e. it adheres 
to the Levy-stable law), and the information about the distribution functions 
if/3, is concentrated in the index a of the stable law. 



3 Probabilistic Interpretation of Empirical Laws 

Since the relaxation rate b cannot be negative, the Levy-stable laws are com- 
pletely asymmetric (supported on the nonnegative half-line) with < a < 1. 
In this case the relaxation function ([6]) with Ajq = N^/°' is well defined and 
takes the KWW form 

/ t \^ 

<i>Kww{t) = Jim (^C{H^r, = exp(- (At)") , (10) 

where A is a positive constant [30|3T] . If « = 1, the relaxation function (11 01) 
becomes (pnit) = exp(— t/r^i) (Debye form), where t^, = is the D relax- 
ation time. Mathematically, this case corresponds to the degenerate case in 
(Q. For any fixed (deterministic) constant A we obtain the only expression 
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( !T0|) . In general, the feature of A is not necessarily true. It is therefore rea- 
sonable to ask what will be, if the constant becomes random. We will get the 
other relaxation laws. 

To find the "scenario" leading to the observable relaxation laws different from 
the KWW form, let us note of the fact that the relaxation function for the 
CC response can write in the form 

oo 

<Pcc{t) = J exp(- {t/Xr) dUX/rcc) = E,{-{t/TccT) , (11) 



where Ea{z) is the one-parameter Mittag-LefHer function, T)j(A/rcc) the one- 
sided Levy-stable probability distribution with the index < a < 1, and tqq 
is constant. It is useful to recall that from the subordination approach the CC 
relaxation response is expressed as 

oo oo 

0cc(t) = / 0D(r)/(t,r)rfr = j exp{-fiT)p^{t,T) dr = E^{-fir) . (12) 



Here the parameter /i is constant, and the contribution of irregular changing 
dipole orientations in the macroscopic evolution of the system is derived from 
the probability density of the directing process S{t). The approach of Weron 
and Jurlewich [2HI25] is based on the other conception. It is that each indi- 
vidual dipole in a complex system relaxes exponentially, but their relaxation 
rates are different and obey a probability distribution (continuous function). 
However, the subordination approach brings advantages in deriving a rather 
simple macroscopic equation for the description of the CC relaxation response. 

The result (fTTl) may be interpreted as a weighted average (or as randomizing 
the parameter A) of the stretched exponential relaxation (fTOj) respect to the 
distribution function Ta(A) of the scale parameter A. This idea works not only 
for the CC relaxation. Really, let Qa be such a random value that its Laplace 
transform is the stretched exponential function 

oo 

(e-'Q-) = J exp{-st)ha{t)dt, < a < 1. (13) 



Then the random value Qa is distributed according to the one-side Levy- 
stable law with the probability distribution function ha{t) with < a < 1 (see 
details, for example, [3lj). Now let the random value Gb be independent of Qa 
and distributed according to the gamma law [32] defined by the probability 
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distribution function 

9bit) = ^t'-'e-\ 6>0, t>0. 



In this connection it should be pointed out that the Laplace transform of Gb 
takes the form 

(e-^ = I exp(-sO gbit) dt = — — ^ . (14) 
^ ^ 

For the random value B Qa (GbY^^ one obtains 

-sBQaGa^^ ^Jexp{-Bst^/'')ha{s)d)j gb{t) dt (15) 

oo 

exp(-(s5)"t)^5(t)rft 



where the positive (arbitrary) constant i? is a scale parameter. The frequency- 
domain response (I)*{uj) is related to the relaxation function 0(t) by the one- 
sided Fourier transform: 

rM=7e-(-^)*. (16) 

^ ^ 



As it is well known [7j, the (dielectric) susceptibility x{^) is directly connected 
with 4'*{u}) by the formula: 



(W) - 



oo 

5 



Xo - Xoo 



where the constant Xoo represents the asymptotic value of x{^)j and Xo is 
the value of the opposite limit. Clearly, the process B Qa {GhY^'' leads to the 
following time-frequency response 

= (i + (Lrr 

This is just the HN relaxation response. Evidently, for 6 = 1 the expression 
corresponds to the Cole-Davidson (CD) empirical law. When the random value 
{GiY^^ follows the WeibuU distribution [32] with the shape parameter equal to 
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c, we arrive at the CC relaxation. Thus, the KWW, CC, CD and HN relaxation 
functions are very close in connection from the probabilistic point of view to 
the random processes associated with the relaxation. 

In this connection it should be pointed out that the evolution of n^{t) and 
rii in Eq. (E]) can be connected with the Mittag-LefHer distribution. Let Z„ 
denote the sum of n independent random values with the Mittag-LefHer distri- 
bution. Then the Laplace transform of is ( l + s"'/n) which tends to 
e"**" as n tends to infinity. Following Pillai [33], this demonstrates an infinity 
divisibility of the Mittag-LefHer distribution. By virtue of the power asymp- 
totic form (long tail) the distribution with parameter a is attracted to the 
stable distribution with exponent a, < a < 1. The property of the Mittag- 
LefHer distribution allows one to develop a corresponding stochastic process. 
The process (called Mittag-LefHer's) arises of subordinating a stable process 
by a directing gamma process [33]. In this case the relaxation function has the 
Havriliak-Negami form 



where a, b, Thn are constant. The one-side Fourier transformation of the re- 
laxation function gives 



oo 



^. (19) 







This result also corresponds to the well-know HN empirical law. Thus, the HN 
relaxation can be explained from the subordination approach, if the hitting 
time process of dipole orientations transforms into the Mittag-LefHer process. 
For that the hitting time process has an appropriate distribution attracted 
to the stable distribution. The subordination of the latter results just in the 
Mittag-LefHer process. It is interesting to observe that the Levy process sub- 
ordinated by another Levy one leads again to the Levy process, but with other 
index [34] . Observe that in this point the subordination approach is almost 
equivalent to the approach studied in |24j. 



4 Internal Structure of Complex Systems 



In any dielectric (complex) system under an week external electric field (ex- 
ternal action) only a part (active dipoles or objects) of the total number N 
of dipoles is directly governed by changes of the field. But even those dipoles. 
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not contributing to the relaxation dynamics, can have an effect on the behav- 
ior of active dipoles. This means that the z-the active dipole interacts with 
Ni — 1 inactive neighbors forming a cluster of size Ni. The number K^^ of 
active dipoles in the system is equal to the number of clusters. The sum of the 
clusters exceeds A^, the size of the system. Because of the screening effects the 
active dipoles can "see" only some of their active neighbors. If so, the coopera- 
tive regions built upon the active dipoles will appear. The number of the such 
mesoscopic regions is determined by their sizes Mi, M2, .... The contribution 
of each region to the total relaxation rate is a sum of the contributions of 
all active dipoles over the region. Generally speaking, the sums are random. 
Hence, the j-th region has its relaxation rate, say (3j]y, equal to 

Mi+...+Mj 
j=Mi + ...+Afj_i + l 



For j = 1 the latter expression is simply the sum Pin = J2iJiPiN- Next for 
j = 2 it takes the form P2N = J^ff^M^i PiN and so on. The relaxation function 
of the whole system 

is provided by the total relaxation rate /^at as the sum of the contributions 
over all cooperative regions: 

Ln 

Pn = ^ (3jN- 



As a rule, the relaxing systems consist of a large number of dipoles so that 
the limit transition /3 = lim7v-s.oo Pn is valid (in practice, N ^ 10^ and more is 
enough). Limit theorems for the random sums have been recently established 
in [35] . 

The number of dipoles directly engaged in the relaxation process is random 
as well as their locations are random too. Obviously, all the quantities A^j's, 
Mj's, /Sjiv's and those defined by them, are random values. Their stochastic 
characteristics determine the total relaxation rate (3^, but they are not known. 
Nevertheless, on the basis of the limit theorems of probability theory, the dis- 
tribution of the limit /3 (for the large relaxing systems) can be defined, even 
with rather information about the distributions of micro/mesoscopic quanti- 
ties. 

In the approach it is quite enough to consider stochastically independent se- 
quences of random values A^i's, M^'s, /Sjat's. Each sequence consists of indepen- 
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dent and identically distributed nonnegative random values that have either 
finite expected value finite or long-tailed distribution. Then the total relax- 
ation rate /3 takes the form corresponding to one of the empirical responses 
(see Section [2]) • It should be noted that the distribution of a nonnegative ran- 
dom value, say X, has a long tail, if and only if the tail Pr(X > x) fulfills the 
condition 

PrfX > x) 

lim — ^ = const > (20) 

x—^oo X~^ 



for some < 7 < 1 so that for value x the tail exhibits the fractional power 
law [17f26]. Many different continuous and discrete distributions are well 
known to satisfy the condition (120]) . Classical examples are completely asym- 
metric Levy-stable laws as well as the Pareto and Burr distributions with an 
appropriate choice of their parameters [17]. To obtain the discrete distribu- 
tions with long tails, one should apply a quantization procedure to the above 
continuous examples [26] . If the distribution of random value X has long tail, 
then the expected value (X) is infinite. The finiteness of the expected value 
and long-tail property ( 120|) can be presented only on different levels (theirs are 
three: an active dipole — )■ a cluster =^ a cooperative region) of the complex 
system. To sum up. Table [1] shows the connection between the internal prop- 
erties of complex system's dynamics and the empirical relaxation responses, 
as well as the physical sense of the parameters characterizing the responses. 
The proposed approach leads to a very general scenario of relaxation, from the 
stochastic nature of microscopic dynamics through the hierarchical structure 
of parallel multi-channel processes to the deterministic macroscopic laws of 
relaxation given by ( ITOj) and f lT7|) . 



5 Energy Criterion 



The common property of the empirical relaxation laws is that they exhibit the 
high-frequency power law in the susceptibility: 

X{uj) oc [ — j for w > cjp. 



where the exponent n falls in range (0,1) and the constant Up is the loss peak 
frequency. As a consequence, for large u the ratio of the imaginary to real 
components of the susceptibility x{^) = x'{^) ~ ix"{^) becomes a constant 
of degree n: 

^ ^ ^ cot (n—\ for u 3> Up. (21) 



X'ioo) V 2 
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Table 1 

The connection between the internal properties of complex systems and their relax- 
ation response (the notations of the column "Parameters" correspond to (jlOp and 
(fT7|) : the constant 7 according to (pO]) ). 



Law 


Parameters 


Ni 






D 


a = 1 
6 = 1 


{Ni) < 00 


{Mj) < 00 


{(3iN) < 00 


KWW 


< a < 1 


{Ni) < 00 


{Mj) < 00 


long tail 

7 = a 


CD 


a = 1 

< 6 < 1 


im) < 00 


long tail 

7 = 6 


{/3iiy) < 00 


CC 


< a < 1 
6 = 1 


long tail 

7 = a 


{Mj) < 00 


long tail 

7 = a 


HN 


< a < 1 
< 6 < 1 


long tail 

7 = a 


long tail 

7 = 6 


long tail 

7 = a 



However, the D response has not the property. It should be noted the physical 
significance of expression (12T]) . At high frequencies the ration of the macro- 
scopic energy lost per radian to the energy stored at the peak is independent 
of frequency. 

Jonscher [2j has advanced a hypothesis that the fact is based on the identical 
property of individual structural elements of the systems. This explains the 
universality in the large scale behavior of complex systems, but needs for the 
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precise derivation. In the framework of the proposed and mentioned-above 
model the physical intuition can be strictly argumentative. Really, the condi- 
tion f l20p applied to any relaxation rate fi leads to the scaling property of the 
relaxation-rate distribution at large h (see also ([H])). The asymptotic behav- 
ior of the distribution is connected with the short-time asymptotic properties 
of the associated relaxation function </'(t), and the response function as its 
derivative f{t) = —d(j){t) / dt takes the form 



for t — )■ 0, where U{t) is a slowly varying function so that limj_^o U{ct)/U{t) = 
1 for any constant c > 0. It may be easily verified that the short-time behavior 
of f{t) corresponds to the high-frequency properties of the susceptibility x{^)- 

XM = x'H - tx"{ou) oc {zco)-^U{l/u). 



The result yields straightforwardly the energy criterion fl21l) with n = 1 — 7. 
The long-tail property of micro/meso/macroscopic relaxation rate with the 
parameter 7 leads to micro / meso / macroscopic energy criterion with the char- 
acteristic constant 1 — 7. The analysis of the model shows [25j that in the 
HN, CC and KWW responses the energy criterion is the case for all mi- 
cro/meso/macro levels, and the constant n for the HN case is defined not 
only the long-tail property of the distribution of cluster sizes, but one of 
cooperative-region sizes. In the CD case the microscopic energy criterion is 
absent. The high-frequency power law of this response results only from the 
long-tail property of the distribution of cooperative-region sizes. 



6 Macroscopic Response Equations 



It is well known (see e. g. [3F]) that the relaxation function 0(t) has to fulfil 
the two-state master equation 

^ = -r(t)0(t), 0(0) = 1, (22) 



where the nonnegative, time-dependent value r{t) is the transition rate of the 
relaxing system (i. e. the probability of transition per unit time). This is a 
macroscopic deterministic equation. Although the equation does not contain 
any (for example, micro/mesoscopic and so on) details about relaxation pro- 
cesses, it is convenient for practical purpose because of its simplicity. So, in 
the case of the KWW relaxation r{t) = aAH""'^. For the D relaxation the 
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equation (122|) has the simplest form r(t) = tq = const. Then, both relax- 
ation function and response function f{t) = — d(f){t) / dt satisfy the same 
equation, and their expressions coincide. However, this equivalence is wrong 
for nonexponential relaxations. Following the old notation ([IS], p. 8), the 
response function is a pulse-response function of the polarization. Any relax- 
ation function or "decay function of the polarization" (by definition, see |l6j) 
tends to 1 for t — )■ 0. But some response functions have a singularity in zero, 
as appears, for example, from the CD response function. In general, the re- 
sponse function f{t) can be a solution of the first-order differential equation 
with variable coefficients: 

m+r{t)l^ = 0. (23) 



dr{t) 
dt 



When the relaxation obeys the CC, CD, HN laws, the equations fl22l) and 
(1231) are not just simple because of a sufficiently complicated expression for 
r{t). It may be attempted to transform fl2^ and fl2^ in a integro-differential 
form simpler than the input ones. Really, the way gives some progress. The 
CC relaxation and response functions can be expressed as a solution of the 
fractional differential equation. 

In this case we have 

(Pccit) = E^i-it/rr), 



and the response function is written in the following series representation 

^ (-l)^(t/r)°(i-^)-^ 

JCCyt - 2^ T^r I , M • 

t^o rr[a{l + k)] 



According to the book of Miller and Ross [37], the one-parameter Mittag- 
Leffler function Ea{—{t/TY) satisfies the identity 

j'-'^[DE^{-{t/rr)] = -^E^i-it/rn 



where D denotes the usual differentiation operator d/dt, and the operator 
J^~" is the Riemann-Liouville fractional integral having the form 

1 r 

rx{t) = —— {t~sy-^x{s)ds. (24) 
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From this it follows that for the CC relaxation the relaxation function fulfils 

J^-'^[D<Pcc{t)] = -^<Pccit), (25) 



and the response function does 

D[j'--fcc{t)] = -^ fecit). (26) 

In this connection it should be pointed out that equation fl25p is expressed in 
terms of the fractional integral of the ordinary derivative, and (12^ is in terms 
of the ordinary derivative of the fractional integral. The equations fl25l) and 
fl26l) are fully equivalent to fl22|) and fl23l) with the transition rate 

rcc(t) = -[E^{-{t/rr)]-'D[E^{-{t/rr)]. 

In fact, this is only another formulation of fl2^ and fl23p . in terms of the Green 
function (see more details, for example, in 



With macroscopic equations for the CD and HN responses the situation is 
more intricate. From the above consideration it is seen that the responses are 
results of the (not simple) integral transformations. This means that probably, 
the CD and HN responses satisfy enough complicated macroscopic equations. 
Therefore, it should be given consideration. 

An interesting idea was suggested in [5]. It proceeds from the fact that the 
ordinary equation of exponential relaxation can be written in the form 

exp(-f]ot)^ exp{not)fY){t) = , 
at 

where is constant. By the direct substitution [5] it is easily verified that 
the CD response function is a solution the following equation 

exp{-not)D''[exp{not)fcBit)] = . (27) 

Here the fractional derivative 

d } 

D'^[x{t)] = l/T{l-u)-j{t-srx{T)dT 



is defined as well as in [5]. Next, Eq. (1271) can be again generalized, namely 
exp(-fiot^)/^1exp(not'^)/ge„(t)] = . (2S 
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Now its solution 

Unit) t"-' exp{-non , 



describes both CD (/i = 1 and < u < 1) and KWW (0 < fi = u < 1) 
response functions. 

In [5] it is proved that the macroscopic equation hke {D^ + 0,)"^^ f(x) = 
is well for the HN response. Their conclusion is based on the proof of the 
operator relation (Appendix A) 

exp{-nuD^-')D'' exp{nuD^-') = {D' + < £ < 1, a < e, 



where m is a variable, and fl a constant (all the notations follow [5]). However, 
the commutator (A2) in [5] does not hold true for any continuous functions. 
In particular, by direct calculations one can found 

[D", ntD^-'] t"-ir(a) = n{a + €-l) -— 

-L \S I 



rather than aQt^~^ /r{e) that it follows from [5]. The cause lies in vanishing 
the term t"~^ = 0. On the other hand, if in time domain one expands the 
Havriliak-Negami (Cole-Cole) response function as an infinite power series, its 
first term is just t°'~^/r(a) (in notation of |5]). Thus the expression (A5) is 
invalid for the HN response. The commutator (A2) remains true only for the 
case, when e = 1 (CD relaxation). The expression (A5) for the Cole-Davidson 
case can be derived directly by the Leibniz's formula for fractional derivatives 
without the operator identity. 

In this connection it should also be recalled about difficulties with the inter- 
pretation of fractional operators in terms of the fractal Cantor set. Following 
[5|39] . the memory function is represented by a Cantor fractal function. How- 
ever, such a memory function possesses only a power-like property asymptot- 
ically, and the approach itself requires else some average procedure for the 
log-periodic term appearing together with the power [10]. The formalism [31] 
is exactly macroscopic, but not stochastic. The main feature of the complex 
system evolution hes just in a stochastic background of dynamical processes. 
In the framework of the approach [39] , to find the relationship between the lo- 
cal random characteristics of complex systems and the universal deterministic 
empirical laws is hardly possible. 
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7 Concluding Renicirks 



In this paper we have shown the outlook of the probabihstic approach proposed 
to the analysis of relaxation phenomena in the complex systems. The approach 
permits ones to consider the observable relaxation law on the unique theoret- 
ical base, the limits theorem of probability theory. The general probabilistic 
formalism treats the relaxation of the complex systems regardless of the pre- 
cise nature of local interactions. In a natural way, it gives an efficient method 
for calculating the dynamical evaluating averages of the relaxation processes. 
We have obtained all (known up to now) the empirical relaxation laws, char- 
acterized their parameters, connected the parameters with local random char- 
acteristics of the relaxation processes, reconstructed the internal structure of 
relaxing systems, justified the energy criterion, demonstrated the transition 
from the analysis of the microscopic random dynamics in the systems to the 
macroscopic deterministic description by integro-differential equations. As a 
rule, the classical methods of statistical physics take into account the Central 
Limit Theorem in respect to the probability distributions having finite invari- 
ance. However, the assumption does not help to clarify the nature of relaxation 
phenomena. The above approach has the advantage over the traditional mod- 
els and goes behind the classical statistical physics. The preliminary results 
presented in this paper are promising and give a confidence for a fundamental 
understanding of the relaxation processes in the framework. 
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